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We study operators acting on a tensor product Hilbert space and investigate their product nu- 
merical range, product numerical radius and separable numerical range. Concrete bounds for the 
product numerical range for Hermitian operators are derived. Product numerical range of a non- 
Hermitian operator forms a subset of the standard numerical range containing the barycenter of 
the spectrum. While the latter set is convex, the product range needs not to be convex nor simply 
connected. The product numerical range of a tensor product is equal to the Minkowski product of 
numerical ranges of individual factors. 



I. INTRODUCTION 



Let X be an operator acting on an iV-dimensional Hilbert space 'Hat. Let A(X) denote its numerical 
range, i.e. the set of all A such that there exists a normalized state ji/;) e Hn, WMl — which satisfies 

In this work we study an analogous notion defined for operators acting on a composite Hilbert space with 
a tensor product structure. Consider first a bi-partite Hilbert space, 

UN^nK^nM, (1) 

of a composite dimension ~ KM . 

Definition 1 (Product numerical range) Let X he an operator acting on the composite Hilbert space 
We define the product numerical range I^{X) of X, with respect to the tensor product structure of 
T-Ln, as 

A®(X) = {(V-A ® ^JsmA <E> ^b) ■■ \4'a) e Hk, \4'b) e Hm} , (2) 
where |Vm) G 'Hk o,nd IV's) G "Hm are normalized. 

Definition 2 (Product numerical radius) Let Hn — Hk ® Hm be a tensor product Hilbert space. We 
define the product numerical radius r'^{X) of X , with respect to this tensor product structure, as 

r'^{X)^ma.x{\z\:zeA^{X)}. (3) 

The notion of numerical range of a given operator, also called "field of values" [U Chapter 1] , has been 
extensively studied during the last few decades [5H3] and its usefulness in quantum theory has been em- 
phasized in]. Several generalizations of numerical range are known - see e.g. [I] Section 1.8]. In particular, 
Marcus introduced the notion of decomposable numerical range [6l[7], the properties of which are a subject 
of considerable interest [8l |9] . 

The product numerical range, which forms the central point of this work, can be considered as a particular 
case of the decomposable numerical range defined for operators acting on a tensor product Hilbert space. 
This notion may also be considered as a numerical range relative to the proper subgroup U (K) x U{M) of 
the full unitary group U {KM) . 

In papers |10fll3| the same object was called local numerical range in view of notation common in quantum 
mechanics. This name seems to be natural for the physicists audience, but to be more consistent with the 
mathematical terminology we will use in this paper the name product numerical range, although a longer 
version "local product numerical range" would be even more accurate. 

In a recent paper of Dirr et al. [TD] some geometric properties of the product numerical range and product 
C-numerical range were investigated. Another paper of the same group [T^] demonstrates the possible 
application of these concepts in the theory of quantum information and the theory of quantum control. 
Product numerical range of unitary operators was very recently used by Duan et al. to tackle the problem of 
local distinguishability of unitary operators [TT]. Knowing product numerical range of a Hermitian operator 
one could solve other important problems in the theory of quantum information, as establishing whether a 
given quantum map is positive, or obtaining bounds for the minimum output entropy of a quantum channel 

The main goal of this paper is to stimulate research on product numerical ranges. We derive several 
bounds for the product numerical range of a Hermitian operator defined on a space with a two-fold tensor 
structure, which corresponds to a bi-partite physical system. In the non-Hermitian case we use the relation 
between the product numerical range of a tensor product and the Minkowski product of numerical ranges 
to establish a general bound for the product numerical range based on the operator Schmidt decomposition. 
We show that a tensor product of two operators, acting on a Hilbert space with a two-fold tensor structure, 
has a simply connected product numerical range. A similar property does not hold for operators acting on 
a space with a larger number of factors. We introduce a class of product diagonalizable operators, for which 
a convenient method to parameterize product numerical range is proposed. 
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Although this work leaves several problems related to product numerical range unsolved, we believe it 
could point out directions for further mathematical research, which will find direct applications in the theory 
of quantum information. In order to invite reader to contribute to this field we conclude the paper with a 
list of exemplary open problems. 

II. PROPERTIES OF PRODUCT NUMERICAL RANGE 

In this section we are going to consider arbitrary operators acting on a bipartite Hilbert space ([T]). 

A. General case 

It is not difficult to establish the basic properties of the product numerical range which are independent 
of the partition of the Hilbert space and of the structure of the operator. We list them below leaving some 
simple items without a proof. 

1. Basic properties 

We begin this section with some simple topological facts concerning product numerical range for general 
operators. 

Property 1 Product numerical range forms a connected set in the complex plane. 

Proof. The above is true because product numerical range is a continuous image of a connected set. ■ 

Property 2 (Subadditivity) Product numerical range is subadditive. For all A, B Cz M„ 

A^{A + B)cA^{A)+A^{B). (4) 
Property 3 (Translation) For all A e M„ and a eC 

A^{A + al)^A^{A) + a. (5) 
Property 4 (Scalar multiplication) For all A e M„ and a G C 

A^{aA) = aA^{A). (6) 
Property 5 (Product unitary invariance) For all A E Mmn 

A^{{U(g,V)A{U(g,V)^) ^A^{A), (7) 

for unitary U £ M,„ and F G M„. 
Property 6 Let A e M„i and B e M„ 

1. If one of them is normal then the numerical range of their tensor product coincides with the convex 
hull of the product numerical range, 

A(A® B) = Co(A^(A(8)B)). (8) 

2. If e^^ A is positive semidefinite for some 9 £ [0, 2ti), then 

A{A®B)^^{A®B). (9) 
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Proof. This property can be proven using Lemma [T] stated in the following subsection and Theorem 4.2.16 
in[T]. ■ 

Let H{A) = 5(A + At) and S{A) = 1{A- At). 
Property 7 (Projection) For all A e M„ 

A^{H{A))^Re A^{A) (10) 

and 

A^{S{A))=ilmA^{A). (11) 
Property 8 The product numerical range does not need to be convex. 
Proof. Consider the following simple example. 
Example 1 Let 



. ,io\ /io\ ./oo\ /oo 
oh +no ih 1 



(12) 



Then A is normal matrix with eigenvalues 0,1, i. It is easy to see that 1 € A** (A) and i G A^{A), but 
(1 + i)/2 ^ A* (A). Actually, by direct computation we have 



A^{A) = {x + yi -.0 < x,0 <y,y/x + y/y <l} . 
Product numerical range of matrix A is presented in Figure [7} 



(13) 




FIG. 1. The comparison of the numerical range (gray triangle) and the product numerical range (dashed set) for 
matrix A defined in Eq. (12 1. 



Product numerical range forms a nonempty set for a general operator. In particular it contains the 
barycenter of the spectrum. 
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Property 9 Product numerical range of A ^km includes the barycenter of the spectrum, 

^ tvA e A® (A). (14) 



KM 

Proof. Let A be an operator acting on a tensor product Hilbert space Hk ® T~Lm ■ Let us write 

1 1 1 1 

iO? =KM^' ^^^^ ® 1)) = ]^ E M^r {A (1 ® |^.)(^.|)) , (15) 

i=l 



where {"iAilf^i is an arbitrary orthonormal basis in "Hm- The last sum in (15 1 is a convex combination of 
elements in the numerical range of tri^, where tri denotes the partial trace with respect to Hk- Remember 
that A (triA) is convex. Hence there exists an element tp S Hm of norm one such that 

M 

J2 ^tr(A(l® = (^ItriAIV) = tr(A(l0 |^)(V'|)). (16) 

i=l 

By repeating the same trick, we can replace the remaining identity with a single one-dimensional projector 
and obtain 

^ tr^= ^tr(yl(l® |V')(V|)) (17) 



KM K 

= tr (A (10) (01 IV) im = (0 ® V'l^l'/' ^ 



for some \4>) G Hk and j^/)) ^Hm of norm one. The last equality in (17) means the same as {tTA)/KM G 
N^{A), which we wanted to prove. ■ 

In the particular case tvA — 0, Property [9] has already been used in [TT]. The above reasoning can be 
generalized to the multipartite case (cf. Section p^II A| . 

Note that the barycenter does not have to lie in the interior of the product numerical range. For example, 
point jj, denoted by the black cross in Figure [ij 

Property 10 Product numerical radius is a vector norm on matrices, but it is not a matrix norm. Product 
numerical radius is invariant with respect to local unitaries, which have the tensor product structure. 

Proof. The only hard part of the proof is positivity. In the Hermitian case this follows from Proposition |3] 
the non-Hermitian case can be reduced to Hermitian by the projection property. I 

Property 11 If A Cz M/^m can be diagonalized to S using product unitary matrices U G M^, V G M^v/ (^i.e. 
there exist unitary U, V such that {U ® V) A (U^ ®V^) = then 

A^{A) ^ {z : z ^ ((pi,P2, • • • ,Pk) fXi (91,^2, • • • ,9m)) • (£1,1,^2,2, • • ■ ,'^km,km)}, (18) 
where J^kPk = 1. J2m 1m 1 and Pk,qrn>0- 

Proof. This follows from more general Proposition |6] in Section |III A[ ■ 

Example 2 We can apply the last property to Example [1} The matrix in this example is given by 

^ = diag(l, 0,0,0- (19) 
By Property \l 1\ we have the following parametrization of the product numerical range of A: 

pq + i{l-p){l-q), p,gG[0,l]. (20) 
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2. Relation to Minkowski geometric algebra 



We shall start this section by recalling the Minkowski geometric algebra of complex sets as developed by 
Farouki et al. 15i. For any sets Zi and Z2 in the complex plane one defines their Minkowski sum, 



Zi S Z2 = {z : z = zi + Z2, zi € Zi, Z2 € Z2} (21) 

and Minkowski product, 

ZiMZ2 = {z:z = Z1Z2, zi e Zi, Z2 e Z2} . (22) 

Note that the above operations are not denoted by ® and (X) as in the original paper p[5], in order to avoid 
the risk of confusion with the direct sum of operators or their tensor product. 

A simple lemma concerning the Minkowski sum and product has interestingly deep consequences. Let us 
define the Kronecker sum of two operators as 

A® B = A®1 + 1(E) B. (23) 

Lemma 1 Product numerical range of the Kronecker product of arbitrary operators is equal to the Minkowski 
product of the numerical ranges of both factors, 

A^{A(»B) = A{A)MA{B), (24) 

while product numerical range of the Kronecker sum of arbitrary operators is equal to the Minkowski sum of 
the numerical ranges of both factors, 

A^{AQ)B) = A{A)mA{B). (25) 
Proof follows directly from the definition of the product numerical range. First part of the above lemma 



has already been used in [TU]. Observe that the definition (22 1 and the property (24) can be naturally 
generalized to an arbitrary number of factors. 

Thus the problem of finding the product numerical range of a tensor product can be analyzed by checking 
what sorts of subsets of the complex plane one can obtain by multiplying two or more numerical ranges. 
This very problem has recently been investigated in a series of papers by Farouki et al. - see [151 1161 and 
references therein. For instance, the structure of the Minkowski product of several discs in the complex plane 
was analyzed in detail by Farouki and Pottmann 16 . 

The above results concerning Minkowski product can be used directly to find the product numerical range 
of a tensor product of an arbitrary number of factors acting on two-dimensional subspaces. The numerical 
range of any matrix of order two forms an ellipse |17j . which may degenerate to an interval or a point. For 
instance it is known that the numerical range of the matrix X = (g '^J'), forms a disk of radius \r\ centered 
at c. Consider a family of operators with a tensor product structure 

Yir^,r2)^X^®X2^(l ) ® ( J ) • (26) 



The product numerical range of Y takes different shapes depending on the values of the radii ri and r2 
of both discs. According to the results of [T^], one may find the values of these parameters for which the 
boundary of the product numerical range is a cardioid, a limagon of Pascal, or the outer loop of a Cartesian 
oval - see Figure [2] 

Analysis of the Minkowski product of two sets becomes easier if none of them contains 0. In such a case one 
may use the log-polar coordinates in the complex plane and reduce the Minkowski product to a Minkowski 
sum in the new coordinates. 

Let us now consider the opposite case. 

Lemma 2 // the numerical range of one of the factors Ai or A2 contains 0, then the product numerical 
range of the tensor product I^{Ai (g) A2) is star-shaped. 
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FIG. 2. Product numerical range for the operator Y defined in (261 with (ri,r2) equal to a) (1, 1) (cardioid), b) 
(0.7, 1) (limagon of Pascal), and c) (0.5, 1.2) (Cartesian oval). 



Proof. We have 



A^iAi®A2)^ y zA{Ai). (27) 

zeA{A2) 



Without loss of generahty, we may assume that g A{Ai). It is known that the numerical range A{Ai) is 
convex. Hence zA{Ai) is star-shaped with respect to for arbitrary z E C We get that A^{Ai (8)^2) is 
star-shaped with respect to and therefore simply connected. I 

It is conceivable that the assertion of Lemma[2]may also hold without the assumption G A {Ai) U A (^2), 
but so far, we were not able to prove or disprove this. 

We are now in a position to formulate the main result of this section. 

Proposition 1 Let Ai , A2 be arbitrary operators. The product numerical range of Ai (g) A2 is simply con- 
nected. 

The above proposition, proved in|^ does not hold for a three-fold tensor product. In this way we confirm 
the conjecture formulated in |10j that one needs to work with at least tripartite systems to construct a tensor 
product operator whose product numerical range is not simply connected. 

3. Inclusion properties 

Generically, operators acting on a bipartite Hilbert space ([T]) do not exhibit the tensor product form. 
However, treating operators as vectors in the Hilbert-Schmidt space endowed with the Hilbert-Schmidt 
scalar product, {A\B) = TiA^B, we may use the operator Schmidt decomposition. In close analogy to 
regular Schmidt decomposition, any operator X acting on Hk T-Lk can be decomposed as a sum of not 



X = VmTAi ®Bi + h ./Jii^iAk'^ ® . (28) 

To find the explicit form of this decomposition, it is convenient to use the reshuffled matrix, Y — X^, such 
that in the product basis it consists of the same entries as the original matrix X, but ordered differently, 
{i,j\Y\i',j') = {i,i'\X\j,j'). Then the non-negative components Hi of the Schmidt vector are equal to the 
singular values of the non-negative matrix YY^ of order K^, while operators Ai and Bi with i = 1, . . . 
are obtained by reshaping eigenvectors of Hermitian matrices, YY^ and Y^Y, respectively [18^. 



Making use of the Schmidt decomposition (28) of an arbitrary bipartite operator X we can formulate 



a proposition concerning its product numerical range. 
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Proposition 2 



A®(X) c ^{A{Ai)MA{Bi))m...m^/jr]^{A{AK2)MA{BK2)) 

= A^iy^Ai (g) Bi) ffl . . . ffl A^{y/fIj^AK2 ®Bk2). (29) 

Proof. It follows directly from Property [2] and from the fact that the product numerical range of a Kronecker 
product of operators equals the Minkowski product of their numerical ranges (Lemma [T]) . ■ 

B. Hermitian case 

1. Baste properties 

Consider a Hermitian operator X = X^ with ordered spectrum Ai < A2 < • • • < Aat. A closed interval 
[Afc, Afc+i] will be called a segment of the spectrum. Thus the spectrum consists of — 1 segments, some 
of which reduce to a single point, if the spectrum is degenerated. The numerical range of a Hermitian 
operator reads A(A") = [Ai,AAr] [I]. In this case also product numerical range is given by an interval, 
A® (X) = [A^jjj, A®jjx]i bounded by the points A^j^^ and A^^x ~ the extremal points of product numerical range 
under the set of all product states. Hence, these points determine the maximal (the minimal) expectation 
values of an observable X among all product pure states. 

Observe that these extremal values determine the product numerical radius, 

r^{X) = max{|A^,„(X)|, |ALx(^)I}- (30) 

This relation simplifies for quantum states. Their positivity implies that r®{p) = A^a^x(p)- 

Note that the product numerical range by definition depends on the concrete form of the tensor product 
Hilbert space 'Htv = 'Hk®'^m- For instance, if X also possesses the similar product structure, X = Xa®Xb, 
then all its eigenstates, — I, . . . ,N, are of the product form and are called separable states. Any pure 

state which is not of the product form is called entangled. Thus in this particular case both ranges are 
equal, I^{X) = A{X) = [Ai, Aat]. This is also the case if the eigenstates and I^at), corresponding to the 
extremal eigenvalues Ai and Aat, are of the product form. Hence the product structure of X is not necessary 
to assure that both ranges coincide. 

Let us consider an arbitrary orthonormal product basis, \i,j) = \i) ® \j) of Hn- The states \i) e Hk, 
i = 1, . . . , K and |j) G "Ha/, « = 1, • ■ ■ , M satisfy the orthogonahty relation = Si^ySjj'. Let us also 

use the composed indices denoted by Greek letters, fi — = to represent X in the product 

basis 

X,, = {^i\X\y) = {^,J\X\^\J'). (31) 

In such a product representation any diagonal element X^^ denotes the expectation value of X in a product 
state l/i) — |j, j), so it belongs to the product numerical range, G M^{X). 

Interestingly, some information about the product numerical range can be obtained even without specifying 
a concrete tensor product structure in Hat and a product basis. 

2. Invariant features 

In this section we investigate some basic properties of the product numerical range of a Hermitian operator, 
X = X\ which hold independently from the splitting of the Hilbert space T-Ln = T^k ®T-Lm- In the general 
case of an operator X acting on a iiT x M space, not so much can be said about its product range. 

Proposition 3 Product numerical range forms, by definition, a non-empty subset of the numerical range, 
7^ A*(X) C A{X). Furthermore, A* (AC) reduces to a single point A if and only if the operator X is 
proportional to the identity, X = AI. Speaking in terms of the standard Lebesgue volume of an interval we 
arrive at the following statement. // Vol[A(X)] > 0, then Vol[A*(AC)] > 0. 
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The assertion that N^{X) is a single point if and only ii X = XI holds also when we do not assume 
Hermiticity of X. 



Proposition 4 For any Hermitian operator X , acting on an N -dimensional Hilhert space, its product nu- 
merical range P^{X) is convex and forms an interval of the real line. 

Proof. Let us assume that oi and a2 belong to N^{X). Hence there exist two pairs of product vectors, 
such that ai = {xi,yi\X\xi,yi) and 02 — (2:2, y2|-^|a;2, 2/2)- Since it is possible to find a continuous family 
of vectors \4>a{s)) £ Hk, which interpolates between |a;i) and \x2) and another family \(I>b{s)) G Hm, 
which interpolates between \yi) and \y2)-, the expectation values of X between these product states interpo- 
late between ai and 02. Thus the entire interval belongs to the product numerical range, [ai, a2] C JSP{X). ■ 

Note that the above reasoning, applied to an arbitrary (non-Hermitian) operator X, shows that the product 
numerical range forms a connected set in the complex plane. However, as shown in further sections of this 
paper, this set does not need to be convex nor simply connected. 

Before we turn to more specific theorems, let us invoke a useful result. 

Lemma 3 Consider a tensor product complex Hilhert space Hpf — Hk ® ■ Then we have the following: 

a) any subspace Sd CZT-Cn of dimension d = {K — 1)(M — 1) + 1 contains at least one separable state, 

b) there exists a subspace of dimension d — 1 = {K — — 1) which does not contain any separable 



Proof. Part a) of the above Lemma follows directly from Proposition 6 in a paper by Cubitt et al. [19J , but 
this fact was earlier proven in [5D1 . 

Concerning part b), for given integers K, M we give an example of family of (A' — 1) x (A/ — 1) matrices 
A^*-*) i = I, . . . , K -I and j = I, . . . , M -I such that no linear combination of {A^*^) : I < i < K,l < j < M) 
is of rank one. 

Let = +£;(i+iJ"+i) for 1 < i < AT and 1 < j < M. Expficitly, let us denote by E^^i'^ the K x M 
matrix containing 1 at the position {i,j) and zeros elsewhere. Suppose that T — 'Yl!i=\ ^Yl!f=\ ciij-^^^^^ is of 
rank one. As the top-right entry of T is zero, then the first row or the Af-th column of T is zero. Thus we 
have ai J = for all j or oli^m = for all i. Deleting the first row or the Af-th column of T we can proceed 
by induction on K and M to get the desired result. 

Thus the subspace spanned by the vectorizations z ^il''^ l^>0 t^^^^ matrices does not contain 

any product state. ■ 

A subspace containing no product states is called completely entangled and it was shown in [22 that 
a generic subspace of dimension (A' — l)(Af — 1) possesses this property. 

In the proof of part b) of Lemma [sj we used the fact that any pure state in a = KM dimensional 
bipartite Hilbert space can be represented in terms its Schmidt decomposition, 



We have assumed here that K < M and denoted a suitably rotated product basis by \i') (g) \i"). The 
eigenvalues fXi of a positive matrix AA^ are called the Schmidt coefficients of the bipartite state The 
normalization condition — (V'lV') — 1 implies that ||^||hs = trAA^^ = 1, so the Schmidt coefficients /i^ 
form a probability vector. 

The state \'ip) is separable iff the AT x M matrix of coefficients A is of rank one, so the corresponding vector 
of the Schmidt coefficients is pure. Hence Lemma [3] is equivalent to the following. 

Lemma 4 Consider a set of k complex rectangular matrices Ai of size AT x M , which are orthogonal with 
respect to the Hilbert- Schmidt scalar product, {Ai\Aj) = TrAlAj = Sij. 



state. 



K M 



K 




(32) 
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a) If k > d = (K — l)(Af — 1) + 1 there exists a complex vector c determining a linear combination of 
these matrices, ylav ~ X]i=i CiAi, such that is of rank one. 

b) Moreover, simple dimension counting arguments presented in JT^ imply that, for a generic set o/(d— 1) 
such matrices Ai, a rank one linear combination does not exist. 

Proposition 5 (General K x M problem) The following lower bound for the product numerical radius 
is true, 

and a symmetric upper bound for the product minimum holds, 

A^i„ < \k-i){m-i)+i- 

Furthermore, there exist Xi = X\ and X2 = X\ acting on Hkm such that A^ax(^i) < ^k+m{Xi) and 

Amin(-'^2) > >^(K-1){AI-1){X2)- 

Before proving this proposition let us state some of its implications. For any Hermitian operator X 
acting on the 2x2 Hilbert space its product numerical range contains the central segment of the spectrum, 
[^2,^3] C A'^{X). A similar statement holds for the 2 x to problem, [Am,Am+i] C A'^{X). Proposition [s] 
implies that for any X acting on a 3 x 3 system A5 S A'^{X). 

Proof. Lemmajsjimplies in this case that any d = (M— l)(_ftr — 1) + 1 dimensional subspace contains a product 
state. Applying this to the subspace spanned by the {M — 1){K — 1) + 1 eigenstates corresponding to the 
(M-l)(i^-l) + l smallest eigenvalues we obtain a bound A^i„ < X(m-i){k-i)+i- The {M - 1){K - 1) + 1- 
dimensional subspace spanned by the eigenstates corresponding to largest eigenvalues also contains at least 
one product state, which implies that A^^x ^ ^k+ai-i- 

To demonstrate the optimality, it is enough to find a single operator X of size KM such that Xi^k-i){m~i) ^ 
A®(X). By Lemma [sj there exists a subspace IC of dimension {K — 1){M — 1) which does not contain 
a separable state. Let • ■ ■ , \'ip{M-i)(K-i)) be some orthonormal basis of this subspace, and let vectors 
\^{M-i){K~i)+i) ^ • ■ ■ I |'0ifAf) enlarge the previous system to the orthonormal basis of "H^Af- We define matrix 
X as 

KM 

X= (33) 

i=(M-l)(K-l) + l 

Each product unit length vector can be written in the above basis as 

(M~l)(K-l) KM 

\x,y)= E ^^1^')- (34) 

i=l i=(A'/-l)(i<'-l) + l 

Because the subspace K. does not contain any separable state, in the above decomposition of a product state 
there exists i £ {{M + 1){K — 1) + 1, ... , KM} such that fit ^ 0. Thus for any product state we have 

{x,y\X\x,y) > = \m-i)(k-i)- (35) 



Note that it is not possible to obtain complete information about the position of the product numerical 
range relying only on unitarily invariant properties. In general, one has to consider a concrete splitting of 
the composite Hilbert space. 
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III. GENERALIZATIONS 



A. Multipartite operators 

1. Defimtion 

It is straightforward to generalize the notions of product numerical range to spaces with a multipartite 
structure. In this section we study Hilbcrt spaces formed by a tensor product of an arbitrary number of 
factors, 

nN = n„,,(s>---<E>n^„ (36) 

with N — mi . . . TOfc. 

For instance, in the definition of product numerical radius one has to perform the maximization over the 
direct product group of local unitary transformations, 

C/ioc := f/(TOi) X U{m2) X • • • U{mk) , (37) 

embedded in the group of global unitaries, U{N) = U{mi . . . rrik). 

2. Parametrization 

The following proposition gives a parametrization of the product numerical range for operators which can 
be diagonalized by a product of unitary transformations. 

Definition 3 (Product diagonalizable operator) We call operator A g Mmim2...mfc product diagonaliz- 
ahle iff there exists unitary operators Umi G , Um^ G , . . . , C/mfc € such that 

{U,m (8) ?7mi «) . . . (8) U^^)A{Urm (8) (8) . . . (8) C/m = S (38) 
and is a diagonal matrix. 

In the special case where the operator A is Hermitian and positive, such states are called classically correlated 
[23j or locally diagonalizable |18j . 

Proposition 6 Let A he a product diagonalizable operator on 'Hmirn2...mk ~ (8 'Hm2 8) ... 8) "Hm^ . By 

^hd2,...dk ~ {hj2, ■ ■ ■ , lk\^\h,l2, ■ ■ ■ Jk) (39) 

we denote the vector of diagonal elements of A after diagonalization using product unitary matrices. With 
the above notation, we have the parametrization of the product numerical range of A, 

{mi — lm2 — 1 m^ — 1 ^ 
- E E ••• EplX^--pirA.,,„...J, (40) 

(r) (r) (r) (r) 

where p}/,... ,_p„^_i > and Pq + ' " + Pmi-i ^ 1 /^^ r ~ 1, . . . ,k. 

Remark 1 Note that if A is a product diagonalizable operator on 'H2'' = ^2 8) ^2 (8 ... 8) 7^2 , then the above 
parametrization simplifies and we have 

K®{A) = {z: z - ({pW,l-p(i)}(8...® • A, 

e [0,1]}, (41) 



where '• ' denotes the scalar product. 

A proof of the above proposition is given in [B] Later we will use this proposition to study some concrete 
examples. 
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3. Properties 



Proposition 7 Product numerical range of an operator X acting on a tensor product Hilhert space H jv 



) 'Hmfc includes the harycenter of the spectrum i.e., 

1 



N 



(42) 



Proof. We can follow the same line of reasoning as in the proof of Property [9j 



The next example, adopted from |12| . shows that, in the case of a three-partite system, product numerical 
range is not necessarily simply connected. 




(a)A® of matrix given by (US 



FIG. 3. Figures obtained from parametrizations (45 1 and (46 1. In both cases the numerical range forms an equilateral 
triangle plotted in grey. 



Example 3 Consider two unitary matrices Ui and U2 written in the standard computational basis {|000), |001), . . . , 

2i7r 2i7r 2i7v 2i7r 2i7r 2i7r 

Ui = diag(l,eT",eT",e~^,eT",e~T",e~T",l) (43) 

and 

U2 = diag(l,eT^,e^,eT^,e"^,e"T^,e"T^,l). (44) 

Both Ui and U2 have identical eigenvalues, however they have different eigenvectors associated with these 
values. 

Product numerical range of the operator Ui acting on a three-partite system is not simply connected, 
as shown in Fig. 3(a)| However, exchanging the position of two middle eigenvalues one obtains a unitary 
operator U2 with the same spectrum, for which the product numerical range is convex and coincides with 



the standard numerical range, I^{U2) = A{U2), see: Fig. 3(b) 

Making use of Proposition [6] we can explicitly parameterize the product numerical range of these matrices. 



'iUi) = [z:z = ({p, l-p}®{q,l~q}(^{r,l- r}) ■ ({C/i}.,)ti} 



(45) 
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and 

A^iU2) = {z:z = ({p, 1 - ri {g, 1 - g} {r, 1 - r}) • i{U2}u)ti} (46) 

for p,q,re [0, 1]. 



E 0.0 






0.0 

Re 



FIG. 4. Product numerical range of matrix A defined in (47 1. Panel a) shows an image of the product numerical range 



(plotted in gray) obtained by parametrization (481, while panel b) shows its image obtained by random sampling. 



The standard numerical range forms a regular octagon. 



Example 4 (Genus 2) Consider a unitary operator of order 16, acting on a four-qubit system. Assume 
that the corresponding matrix A is diagonal in the product basis {|0000), . . . , |1111)}, 

A = diag(eT ^i^i^e^ ,—l,e~^ ,e^^ ,1,-1, ,e~^ ,l,e^ , —i, — z, e^). (47) 
Applying Proposition [6j we parameterize the product numerical range of this matrix in the following manner 
k®{A) = {z : z = ({p, l-p}®{q,l-q)®{r,l-r}® {s, I - s}) ■ a} (48) 

for q, r, s G [0, 1], where '•' denotes the scalar product. Product numerical range of the matrix A acting on 
the four-partite system forms a set of genus 2, as demonstrated in Figure |4] 

IV. CONCLUDING REMARKS 



In this work we established basic properties of product numerical range (PNR) of operators acting on the 
Hilbert space with a tensor product structure. Even though the definition of this extension of the standard 
numerical range is rather straightforward, characterizing this set for a general operator occurs to be a difficult 
problem. We have not managed therefore to present its general solution, but we obtained concrete bounds 
for this set (from inside and from outside) useful in various cases. Several special cases of this general 
problem are relevant in view of practical applications of product numerical range in the theory of quantum 
information [T3]. To stimulate interest of the mathematical community in this subject we conclude the work 
by providing a list of open problems and sketching their motivation stemming from physics. 
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1. For a given Hermitian operator H acting on a bi-partite Hilbcrt space, Hn = 'Hk ® find its local 
numerical range [A®i„,Aga^]. 

Also precise bounds for PNR in this case would be useful. Of particular importance would be methods of 
establishing whether for a given H the inequality A^j^^ > is satisfied. This statement is equivalent of saying 
that the operator H is block positive (with respect to the given decomposition of the space "Hat), hence the 
corresponding quantum map is positive [Ml ■ 

2. Find the local product range A'^{U) for a possibly large class of unitary operators acting on a Hilbert 
space with the tensor product structure. 

It is equally important to establish bounds applicable for PNR of an arbitrary unitary U, suitable to answer 
the question if S A'^{U). This very issue occurs to be crucial in solving physical problems related to local 
distinguishability of unitary quantum gates [TTl [2] . 

3. For a given operator X acting on a Hilbert space with the tensor product structure find bounds for 
minimum of the modulus \z\, such that z € A^{X). 

This question is related to finding the optimal gate fidelity [IB] or optimizing local fidelity between two 
arbitrary quantum states |14j . 

4. Analyze the genus of PNR in the general case of operators acting on an fc-partite Hilbert space, 
"Hi (g) • • • (gj'Hfc. 

In particular verify the conjecture that the genus g is not larger than k — 2. If this conjecture does not 
hold in general one may still try to check whether it is true for operators that are fc-fold tensor products. 
This issue is then related to the question what the genus of the Minkowski product of k convex sets on the 
complex plane is. 

As a more specific question, one may ask the following: 

5. Is A® {A (E) B) star-shaped for arbitrary operators A, B acting on T-Lk, T~Lm, respectively? 
We have already mentioned that problem in Section p^I A 2[ 

Furthermore, one may also consider similar problems in a more general set-up by studying corresponding 
product analogues of other generalizations of the notion of numerical range. For instance, mathematical 
results on 'local (product) C-numerical range' and 'local (product) C-numerical radius' [inillllllS], separable 
numerical range |14j and 'Liouville numerical range' [26j will find direct applications in several problems in 
theoretical physics. 
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Appendix A: Proof of Proposition [T] 

As a consequence of Lemma |2j we can limit our proof to the case ^ A {Ai) U A (^2)- Rotations of A {Ai) 
or A(A2) in the complex plane cannot affect the genus of A (^1 ® A2) = A(^i) M K{A2). Thus we may 
assume for the beginning that n (A (^i) U A {A2)) = 0. We can achieve this by suitably rotating A (^i) 
and A (A2), because M+ n A {Aj) = or — K+ n A {Aj) = for either j, as a consequence of the convexity of 
Aj and the assumption that ^ Aj. Now we can rotate A (^1) and A {A2) clockwise, so that = becomes 
the minimal number in [0, 2tt) for which M+e*"^ n A {Ai) 7^ 0, as well as the minimal number for which 
K+e** n A (^2) 7^ 0. By our construction, keeping in mind that the sets Aj are closed, there exist ei, £2 > 
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such that K+e n A {Aj) = 0Vjg(o,ej) for j — 1,2. Using the convexity of A (Aj) for j = 1,2 again, we see 
that both A {Aj) must be contained in the closed upper half-plane. Given K+ n A (Aj) ^ 0, by the same 



argument as above we get 



n A(Aj-) 



' n A {Aj 



i. Closeness of A^ implies that 



max {</. e [0, 2tt) \R+e"^ f] A {Aj) 7^ 0} < tt. 



(Al) 



for j — 1, 2. Since ^ A(Ai) U A{A2) and < 27r, we can describe the sets Aj in the log-polar coordinates 
S (z) = (log |z|, Arg {z)). Because the sets Aj are convex, M+e** n A {Aj) must be a closed nonempty interval 
or a point for arbitrary (p and j. Thus in the new coordinates we can write A {Aj) as 



E{A{A,))= y [ri(0,r;(0] xU} 



(A2) 



Using the convexity and closeness of Aj, it is easy to prove that the functions r^. are continuous, which 
implies that (A2) is a simply connected subset of K^. Because the multiplication of complex numbers 
zi = |zi|e*'^i,Z2 = |z2|e"^^ corresponds to the addition of their arguments 0i,^2 as well as the addition of 
the logarithms of their modules, log|zi| and log|z2|, it is straightforward to write an explicit formula for 
A {Al) A {A2) in the coordinates S, 



u 



u 



K iv) rl {i-^), rl {yp) rl - ^)] x {£} 



or 



U [i?+(o,i?-(0]xU}, 



where the functions arc defined in the following way. 



R+ (0 = maxj^^^(o,c-0^),min(0;,e)] (^+ (<^) ^+ " ^)) 
R- (0 = minj„^^^(o^^_^^)_,^^;„(^^^^)j (ri (^) r2 - (^)) 



(A3) 



(A4) 



(A5) 
(A6) 



It turns out that i?+ and i?_ are continuous and the interval (^) , i?+ (^)] is nonempty for arbitrary 
C S [0, (/)^ + . Therefore (A3) is a simply connected subset of M x [O, cf)\ + . Since + 4>\ < 27r, 
S is a homomorphism between M x [O, (f)\ -I- 0^ and 5"^ (M x [O, (t)\ + (t>\])- Therefore A (^i) Kl A {A2) is 
simply connected, as a preimage of a simply connected set by the isomorphism S. This is what we wanted 
to prove. 



Appendix B: Proof of Proposition [S] 



Each product vector \x) = \xi) (g) \x2) 



where 



Thus we have 



{x\A\x) 



\xk) of unit norm can be rewritten in its computational basis. 



(xi^o|0) H h xi^„ii-i\nii - 1)) (g> . . . (g> (a::fc,o|0) H h Xk,,nk-i\'mk - 1)), 



Fr.ol 



\xr.m. -iP = 1 for r = 1 . . . fc. 



(Bl) 
(B2) 



mi— lm2 — 1 TTifc — Imi — 1 m2 — 1 — 1 

E E--E E E-E 

;i=o i2=a ik=o si=o s2=o «fc=o 

^th^th ■ ■ ■ ^fc.ifc 2:1,^1X2,^2 • ■ • Xk,sk {hh ■ • • ^fc|-4|siS2 • ■ • Sfe). 



(B3) 
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Now we must note that, since A is diagonal with respect to the product computational basis, we have 

{hh ■ ■ ■ lk\A\siS2 . . . Sfc) = Xhh...lJ{hl2...h,s,S2...s^}- (B4) 

Thus we get 

mi — 1 m2 — 1 — 1 

{x\A\x)=J2 E ••• E l^i.iin2;2,/.p...|xfc.,j2A,i,,..,,. (B5) 

;i=o i2=Q ik=o 

If we denote p[^^ :— \xr.i\'^, then we have pl[^ + p''i^ + ■ ■ ■ + pin^-i = 1 and 

mi— lm2 — 1 mfc — 1 

{x\A\x}^Y. E ■ • • E pI'M? • ■ -Plf A,,,,...,. (B6) 

;i=o (2=0 /fe=o 

Now it is easy to notice that if we take all possible product states, we will obtain all possible decompositions 

. . . ,p^ > 0, +pW + . . ■+pI:1 = 1 (B7) 
for r = 1, . . . k. Thus the proof is complete. ■ 
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